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Astratto 


Sebbene gli studi clinici per le applicazioni mediche del dimetilsolfossido (DMSO) 
abbiano riportato tossicità negli anni '60, in seguito la FDA ha classificato il DMSO 
nella categoria dei solventi piu sicuri. II DMSO è diventato ampiamente utilizzato in 
molti campi biomedici e gli effetti biologici sono stati trascurati. Nel frattempo, la 
scienza biomedica si è evoluta verso tecniche sensibili ad alto rendimento e nuove 
aree di ricerca, tra cui epigenomica e microRNA. Considerando il suo ampio utilizzo, 
soprattutto per la crioconservazione e in vitrotest, abbiamo valutato l'effetto biologico 
del DMSO utilizzando queste innovazioni tecnologiche. Abbiamo esposto microtessuti 
cardiaci ed epatici 3D su un terreno con o senza 0,1% DMSO e analizzato i profili di 
trascrittoma, proteoma e metilazione del DNA. In entrambi i tipi di tessuto, l'analisi del 
trascrittoma ha rilevato > 2000 geni espressi in modo differenziale che influenzano 
processi biologici simili, indicando così azioni incrociate coerenti del DMSO. Inoltre, 


l'analisi del microRNA ha rivelato deregolazioni su larga scala dei microRNA cardiaci 
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ed effetti più piccoli, sebbene ancora massicci, nei microtessuti epatici. Anche i 
modelli di metilazione dell'intero genoma hanno rivelato la specificità del tessuto. 
Mentre i microtessuti epatici hanno dimostrato cambiamenti non significativi, | 
risultati dei microtessuti cardiaci hanno suggerito l'interruzione dei meccanismi di 
metilazione del DNA che portano a cambiamenti a livello del genoma. | cambiamenti 
estremi nei microRNA e le alterazioni nel panorama epigenetico indicano che il DMSO 
non é inerte. Il suo utilizzo dovrebbe essere riconsiderato, soprattutto per la 
crioconservazione di embrioni e ovociti, poiché può avere un impatto sullo sviluppo 


embrionale. 


[mme n 


Introduction 


Dimethyl sulfoxide (DMSO) is an organic polar aprotic molecule with an amphipathic 
nature that is ideal for dissolving poorly soluble polar and non-polar molecules. 
DMSO is widely used as solvent in toxicology and pharmacology, for cryopreservation 
of cells, and as penetration enhancer during topological treatments. The use of DMSO 
is so obvious that applied concentrations are often unreported. DMSO is generally 
accepted as nontoxic below 10% (v/v) and, in practice, it is assumed that effects of 
DMSO are negligible'??. 


Research from the 1960's till 1990's pinpointed a lot of biologically relevant effects of 
DMSO. It was even actively investigated for medical use, but because of adverse 
effects especially affecting the eyes, most clinical trials were halted in 1965 by the 
United States Food and Drug Administration (FDA). At this point in time, DMSO was 
seen as extremely toxic, comparable to thalidomide^?. The effects of DMSO differ 
depending on dose and route of administration. Introducing DMSO concentrations 
higher than 5096 into the blood resulted in instant hemolysis, white cell stacking and 
fibrinogen precipitation while direct injection of DMSO intravenously could cause 


local irritation and necrosis?. LD50 values gained from monkeys indicated that 880 
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grams applied on the skin or 320 grams injected intravenously would result in 50% 
mortality in 80 kg humans’. This relatively high dose changed the view on DMSO 
toxicity and the FDA classified DMSO in the same class as ethanol, namely class 3 
solvent, which is the safest category with low toxic potential at levels normally 


accepted in pharmaceuticals?. This made the wide use of DMSO possible. 


While generally used at relatively low concentrations, DMSO still has medically useful 
properties such as inducing anti-inflammation, nerve blockage (analgesia), diuretics, 
vasodilation and muscle relaxation’. Furthermore, in cell biology, DMSO is also used 
as inducer of cell differentiation, free radical scavenger and radioprotectant, but most 
often for cryopreservation. Cell cultures for research are often stored in liquid nitrogen 
using slow cooling methodology. To prevent damage by intracellular ice crystals, cells 
are slowly cooled to -80 degrees in the presence of 10% DMSO before storage in 
liquid nitrogen?. However, this procedure is insufficient for biomedical applications. 
Because slow cooling still induces damage due to extracellular ice crystals, 
cryopreservation of human oocytes and embryos for in vitro fertilization (IVF) is 
accomplished by vitrification, in which higher concentrations of cryoprotectants are 
used to prevent ice formation not only in the cells, but in the entire solution. Higher 
concentrations are achieved by using mixtures (for example 1596 DMSO in 
combination with 1596 ethylene glycol) which reduces the amount and toxicity of the 
individual cryoprotectants ?. Furthermore, DMSO readily crosses most tissue 
membranes of lower animals and man’. For example, 2 hours after topological 
application, DMSO had penetrated all investigated hard and soft tissues (6 and 15 
resp.) in rats! '. DMSO is also able to enhance the permeability of other low molecular 
weight compounds, thereby making them pass membranes or going deeper into a 
tissue as they normally would'?, a property highly useful in topological therapies. In 
2009, the first to obtain FDA approval for topological DMSO usage was PENNSAID®, 
which contains diclofenac in a carrier with 45.5% DMSO'*''^, This relatively high 
concentration of DMSO in topological applications is necessary because the skin is 


harder to penetrate than cell membranes. 


Most insights into the molecular effects of DMSO were obtained last century often 
using high doses. In the meantime, biomedical science has evolved towards more 
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sensitive high-throughput techniques and towards new areas of research, including 
epigenome modifications and microRNA-mediated gene silencing. Considering its 
wide use in many biological fields, we analyzed a relatively low dose of DMSO (0.1%, 
which is commonly applied in cell assays) to study the impact on the proteome, 
transcriptome and the epigenome. For this, we exposed in vitro 3D microtissues (a 
maturing iPSC-derived cardiac model and a mature hepatic model) to 0.1% DMSO 
and we collected samples in triplicate at 7 different time points during 2 weeks 
exposure (2h, 8h, 24h, 72h, 168 h, 240 h and 336 h). Thereafter, DMSO effects were 
assessed with full transcriptome analysis (using ribo-depleted total RNA sequencing 
and microRNA sequencing), whole-genome methylation profiling (using MeDIP-seq) 
and proteomics analysis (using mass spectrometry). Our analysis clearly demonstrated 
that DMSO cannot be considered biologically inert but induces large alterations in 
microRNAs (miRNA) and epigenetic landscape, especially in the maturing cardiac 


model. 


Results 


| microtessuti umani 3D (MT) di un modello cardiaco in fase di maturazione e un 
modello epatico maturo sono stati esposti a terreno di coltura con o senza 0,1% 
DMSO per due settimane con punti temporali di campionamento a 2, 8, 72, 168, 240 e 
336 ore. Il proteoma (circa 2.000 proteine misurate), il trascrittoma completo 
(compresi i MIRNA) e la metilazione dell'intero genoma sono stati misurati su 
materiale ottenuto dallo stesso campione. Figura 1contiene una panoramica grafica 
del disegno sperimentale. Per ottenere una prima panoramica degli effetti cross-omici 
indotti da DMSO, per ciascuna piattaforma vengono riepilogate le quantità di entità 
modificate in modo differenziale (tutte corrette per test multipli utilizzando FDR 
«0,05). Il numero di entità variamente modificate differiva tra i tipi di tessuto, con 
campioni cardiaci che mostravano un effetto maggiore del DMSO rispetto a quello 
epatico, ad eccezione degli mRNA. Questa differenza é particolarmente evidente per i 
miRNA e la metilazione del genoma. Poiché i dati proteomici erano meno informativi 


a causa della loro natura parziale, questi risultati sono stati inclusi nei dati 
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supplementari. Inoltre, l'analisi delle componenti principali (PCA), utilizzando medie di 
triplicati, per ciascuna piattaforma (Fig. 2& Dati supplementari) illustra chiare 
differenze tra lo 0,1% di DMSO esposto (DMSO) e i campioni non trattati (UNTR), ad 


eccezione della metilazione nei MT epatici. 


Figura 1 
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Panoramica grafica del disegno sperimentale combinato con sintesi delle entita differenziali 
di ciascun metodo di analisi. Le informazioni specifiche del tessuto sono rappresentate in 
arancione per il cuore e in verde per l'epatico. Inoltre, le esposizioni sono colorate in blu e le 
piattaforme di misurazione in viola. Abbreviazioni: h = ore; MRNA = RNA messaggero; 


miRNA = microRNA. 


figura 2 


https://www.nature.com/articles/s41598-019-40660-0 


5/39 


24/03/22, 06:09 








a) PCA Cardiac RNA Samples 
11 A 
i |] A 
ilo 
E 
8 e 
£14 
n o 
di ee 
È | A 
$1 
sl 
q a 


PC1: 86% variance 














b) PCA Cardiac miRNA Samples 
8 
a) 
o 
î ae a 
3 
£ 1 ie oe 
Aa 
hd 
z 
4 
1. 
è a. 
$ a 





sen de4^ jx delet de 


PC1:90% variance 


c) PCA cardiac DNA methylation 

è A 

o 
o 

N 
g 
8 o 
$ 
È ©) A 
E 9 
S [e E 

~ a 

" aa 

4 2 0 2 


PCI: 17.91% variance 





PC2: 14% variance 


PC2- 8.9% variance 


PC2: 8.89% variance 





DMSO induce cambiamenti drastici nei processi cellulari umani e nel panorama epigenetico in vitro | Rapporti scientifici 


PCA Hepatic RNA Samples 
« [E 4 DMSO 
, oe e UNTR 
* 024 = 002 
e [|t * 008 
e| * 240 = 024 
— = 072 
^| m 168 
a 4 = 240 
= 336 
e 9 
o 
A à á 
A 





0000€ 300000 D 0000€ tooo 


PC1:80% variance 














PCA Hepatic miRNA Samples 
A 
A 
o 
A 
A 
A 
5 
o 
o 
o 
o 





PCA hepatic DNA methylation 
e 
E ha e ^ 
A 9 
A A 
A 
o 
A 
A 
o 
4 2 0 2 4 


PCI: 10.69% variance 


PCAs depicting differences between DMSO and UNTR for all measured platforms. (a) PCA of 


RNAs indicates clear differences in RNA expression between DMSO (triangle) and UNTR 


(circles). Cardiac samples (left) are more distinct from UNTR than hepatic samples (right). (b) 


PCA of miRNAs reveals clear separation between DMSO and UNTR in cardiac samples, while 


hepatic samples seem more susceptible to the duration of the exposure (as seen by colour 
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pattern that corresponds to the specific time points, see legend). (c) PCA of promotor 
methylation indicates differences between DMSO and UNTR for cardiac samples but not for 


hepatic samples. 


To study the molecular effects of DMSO, affected cellular processes were analysed 
using the full transcriptome. Thereafter, tissue-specific effects of DMSO on regulation 
of gene expression were investigated by analysing changes in miRNAs and genome- 


methylation. 


DMSO effects on cellular processes 

DMSO effects on mRNAs were depicted by of PCA (Fig. 2a,b). The clear separation 
between UNTR and 0.1% DMSO indicated that DMSO was able to affect cellular 
processes by altering gene expression. Comparison between DMSO and UNTR 
resulted in 2051 differentially expressed genes (DEGs = FDR <0.05; of which 871 with 
llog2FC| >1) in cardiac MTs and 2711 DEGs (of which 1879 with |log2FC| » 1) in hepatic 
MTs, of which 60.7% and 62.9% DEGs were downregulated respectively. 


Pathway analysis of DEGs 

Per identificare i processi cellulari interessati dall'esposizione a DMSO, i DEG sono stati 
utilizzati per l'analisi della sovrarappresentazione del percorso utilizzando 
ConsensusPathDB !? con il database Reactome curato 16 . Percorsi significativamente 
sovrarappresentati (valore q < 0, 05) sono stati ordinati utilizzando le connessioni 
gerarchiche tra (sotto) percorsi ottenuti dal browser del percorso del reactome 
(tabelle supplementari 1,2 ). Un riepilogo contenente i livelli di percorso gerarchico 
più elevati (d'ora in poi indicati come cluster) è incluso nella Tabella 1 per i MT 


cardiaci e nella Tabella 2 per i MT epatici. 


Tabella 1 Analisi del percorso di DEG e DEP cardiaci rilevati dopo 
un'esposizione allo 0,1% di DMSO. 
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Tabella 2 Analisi del percorso di DEG e DEP epatici rilevati dopo 
un'esposizione allo 0,1% di DMSO. 


Attraverso l'analisi del percorso sui DEG, sono stati trovati 225 percorsi 
significativamente sovrarappresentati (valore q < 0,05) nei MT cardiaci, che 
corrispondevano a 19 cluster (su un totale di 25 cluster nel Pathway Browser) e 167 
percorsi corrispondenti a 16 cluster in MT epatici. C'era una sostanziale 
sovrapposizione tra i tipi di tessuto, con 60 percorsi e 15 cluster trovati in entrambi. 
Sebbene vi fossero differenze nell'entità dell'effetto del DMSO tra i tipi di tessuto, i 


processi biologici interessati dal DMSO non sembrano essere specifici del tessuto. 


Il cluster più colpito in modo significativo da DMSO è il cluster "metabolismo" nei MT 
epatici. Qui, la maggior parte degli effetti sono stati trovati nei percorsi "ciclo 
dell'acido citrico e trasporto di elettroni respiratori" (valore q: 3,5 * 10-19, 63 DEG su 
171 geni, 76,2% sottoregolato), "metabolismo del glucosio" (valore q: 9,9 * 10 ? , 36 
DEG su 77 geni, 80,5% sottoregolato) e "metabolismo di lipidi e lipoproteine" (valore 
q: 2,9 * 10 6, 165 DEG su 728 geni, 55,2% sottoregolato). | cambiamenti nel 
metabolismo dei lipidi e delle lipoproteine non sono stati rilevati nei MT cardiaci 
utilizzando i DEG, ma sono stati osservati effetti simili del DMSO per "ciclo dell'acido 
citrico e trasporto degli elettroni respiratori" (valore q: 1,3 * 10 -12, 58 DEG su 171 
geni, 65,5% sottoregolato) e "metabolismo del glucosio" (valore q: 2,8 * 10 ? , 20 DEG 
su 77 geni, 55,096 sottoregolato), sebbene meno geni fossero sottoregolati nei MT 


cardiaci rispetto a MT epatici. 


Un altro cluster altamente colpito dal trattamento con DMSO in entrambi i tipi di 
tessuto é stato il "trasporto mediato dalle vescicole". Gli effetti del DMSO in questo 
cluster sono stati rilevati principalmente nei processi relativi al trasporto e alla 
secrezione di proteine mediate da Golgi. Di questo processo, il "trasporto anterogrado 
da ER a Golgi" è stato il più colpito (cardiaco: valore q: 2,4 * 10 © , 38 DEG su 134 geni, 
68,4% sottoregolato; epatico: valore q: 1,0 * 10 ? , 44 DEG su 134 geni, il 44,4% 
sottoregolato). Questo percorso faceva anche parte del cluster "metabolismo delle 


proteine", che ha inoltre rivelato gli effetti del DMSO sulla "glicosilazione N-linked 
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dell'asparagina" (Cardiaco: valore q: 2,1 * 10 © , 64 DEG su 283 geni, 76,6% 
sottoregolato ; Epatico: valore q: 1,2 * 10 7, 83 DEG su 283 geni, il 54,896 
sottoregolato), che erano modificazioni proteiche post-trascrizionali necessarie per il 


trasporto di proteine dall'ER al Golgi !”. 


Sebbene il cluster di "risposte cellulari allo stress" sia stato rilevato in entrambi i tipi di 
tessuto, i percorsi sovrarappresentati differiscono. La "senescenza cellulare" è stata 
significativamente influenzata nei MT cardiaci (valore q: 2,6 * 10 ^ dimensione del 
percorso: 192, 42 DEG, 64,3% sottoregolato), ma non rilevata durante l'analisi del 
percorso nei MT epatici. Infine, ulteriori cluster specifici cardiaci affetti da DMSO erano 
"ciclo cellulare", "riparazione del DNA", "biogenesi e mantenimento degli organelli", 
ma anche il cluster altamente significativo (q <0,01) di "organizzazione della 
cromatina". Poiché i percorsi più colpiti sono stati trovati in entrambi i tipi di tessuto, 


ciò indica azioni robuste del DMSO. 


Effetti del DMSO osservati nella regolazione dell'espressione e della traduzione genica 
| percorsi correlati alla regolazione dell'espressione e della traduzione genica sono 
stati già rilevati durante l'analisi del percorso dei MT cardiaci (Tabella 3 ). Sebbene 
questi percorsi non fossero sovrarappresentati nei MT epatici, le quantità di DEG 
rilevate indicavano che il DMSO era in grado di influenzare questi processi e poteva 


indurre alterazioni biologiche al modello cellulare. 


Tabella 3 Percorsi relativi alla regolazione trascrizionale. 


Effetti DMSO osservati nella regolazione dai microRNA 


L'influenza tissutale del DMSO è stata osservata nei dati di sequenziamento dei 
miRNA. Su 1.105 miRNA cardiaci sequenziati, 704 (= 63,7%) erano espressi in modo 
differenziale (DE, FDR <0,05) con il 59,5% che mostrava una downregulation. Inoltre, il 
grafico PCA (Fig. 2c ) ha rivelato una chiara differenza tra i miRNA dei gruppi di 
trattamento. Nelle MT epatiche, su 1.033 miRNA sequenziati, 186 (= 18%) erano DE 
con circa la meta dei miRNA sovraregolati (47,3%). Inoltre, il grafico PCA non solo ha 


rivelato una chiara separazione tra i miRNA dei gruppi di trattamento, ma ha indicato 
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anche un effetto dipendente dal tempo. Le due componenti principali (tempo e dose) 
rappresentavano insieme piu del 95% della variazione totale, lasciando solo un effetto 


minore per altri fattori. 


Per studiare la fonte della differenza tessuto-specifica, i cambiamenti dell'espressione 
genica nel processo di biogenesi del miRNA sono stati studiati in modo più 
dettagliato (Fig. 3 ). La biogenesi del miRNA inizia con la trascrizione del miRNA 
primario trascritto dalla polimerasi Il (un complesso di 11 subunità) o dalla polimerasi 
Ill (un complesso di 10 subunità). Nelle MT cardiache, la polimerasi Il conteneva 1 geni 
sovraregolati e 8 sottoregolati. Nelle MT epatiche, 5 geni erano sottoregolati e 1 
sovraregolato. Sebbene un minor numero di geni fosse sottoregolato nei MT epatici, i 
cambiamenti di piega erano maggiori rispetto ai MT cardiaci. La scissione del pri- 
miRNA in pre-miRNA non sembrava influenzata dall'esposizione al DMSO. D/ICERT(che 
scinde i pre-miRNA) è stato sottoregolato nei MT cardiaci e sovraregolato nei MT 
epatici, descrivendo una chiara differenza specifica del tessuto nella risposta 
all'esposizione al DMSO. Infine, AGO2 (che codifica per il componente principale del 
MIRNA-RISC) è stato sottoregolato nei MT cardiaci. | cambiamenti in questo processo 
potrebbero indurre differenze nel contenuto di miRNA delle cellule e quindi 


influenzare la loro funzione di regolazione. 


Figura 3 
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Effetto DMSO nel processo di silenziamento genico da parte degli RNA. DEG nella biogenesi 
del miRNA. Il processo completo è suddiviso in sottoprocessi (ovali viola) e raffiguranti geni 
coinvolti (rettangoli blu) e DEG rilevati per campioni cardiaci (rettangoli arancioni) ed epatici 


(rettangoli verdi). 
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Nelle MT cardiache, i geni sottoregolati nella biogenesi dei miRNA potrebbero 
spiegare la quantità estremamente elevata di miRNA DE sottoregolato e indicare 
un'estrema deregolazione del silenziamento genico da parte dei miRNA. Nelle MT 
epatiche, dove l'unico cambiamento significativo era la sovraregolazione di DICER7 


(Tabella 3 ), la biogenesi del miRNA non è stata interrotta dall'esposizione al DMSO. 


Per ottenere un'indicazione degli effetti del DMSO sul silenziamento genico da parte 
dei miRNA, il database miRTarBase 1? , contenente interazioni microRNA-target (MTI) 
convalidate sperimentalmente, è stato utilizzato per ottenere bersagli genici dei 
MIRNA DE rilevati. Abbiamo incluso solo MTI con prove evidenti (convalidate da 
reporter assay, Western blot o qPCR). Dei 704 miRNA cardiaci DE, solo 281 (= 40%) 
sono stati trovati nel database, risultando in un totale di 2051 bersagli genici 
potenzialmente interessati dai cambiamenti indotti da DMSO nel silenziamento 
genico del miRNA. Per le MT epatiche, sono stati ottenuti target per 106 DE miRNA (= 
57%), con un totale di 545 geni potenzialmente interessati. | bersagli genici ottenuti 
sono stati utilizzati per visualizzare i percorsi sovrarappresentati utilizzando il browser 
del percorso del reactome (Figura 3 supplementare). Inaspettatamente, i percorsi 
sovrarappresentati si trovavano negli stessi cluster per entrambi i tipi di tessuto. La 
maggior parte degli effetti è stata osservata per "Trasduzione del segnale", "Sistema 
immunitario" ed "Espressione genica". Tuttavia, l'estrema deregolamentazione degli 
MT cardiaci e le informazioni limitate sugli MTI stanno rendendo irrilevante qualsiasi 


analisi a valle sull'mRNA presunto interessato. 


Effetti DMSO osservati in epigenetica 

Al fine di valutare le alterazioni epigenetiche introdotte da DMSO, ci concentriamo 
sulla metilazione del DNA dell'intero genoma. L'analisi del percorso dell'espressione 
genica cardiaca ha rivelato la deregolamentazione dei percorsi di metilazione del 
DNA. Le metiltransferasi DNMTT , un fattore chiave per il mantenimento della 
metilazione del DNA, e DNMT3A, che facilitano sia il de novo che il mantenimento 
della metilazione del DNA, sono state sovraregolate mentre TET7 , che svolge un ruolo 
chiave nella demetilazione attiva, è stata sottoregolata nelle MT cardiache (Fig. 4a). La 
sovraregolazione degli scrittori epigenetici e la sottoregolazione delle gomme dopo il 


trattamento con DMSO hanno indicato l'ipermetilazione genomica, che 
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potenzialmente ha ridotto l'attività trascrizionale. Al contrario, nell'analisi della via 
epatica, non è stata osservata la deregolazione della metilazione del DNA. Si noti che, 
in entrambi i tipi di tessuto, sono state osservate prove trascrizionali per la 
deregolazione di altri meccanismi epigenetici correlati, come la metilazione dell'istone, 
in cui 16 geni e 13 geni erano espressi in modo differenziale rispettivamente nei MT 


cardiaci ed epatici. 


Figura 4 
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Regolazione epigenetica dell'espressione genica. ( a ) DEG coinvolti nella metilazione del 
DNA. Il processo é suddiviso in sottoprocessi (ovali viola) e raffiguranti geni coinvolti 
(rettangoli blu) e DEG rilevati per campioni cardiaci (rettangoli arancioni) ed epatici 
(rettangoli verdi). ( b ) Arricchimento relativo dei DMR (DMSO vs UNTR, FDR 5%) all'interno 
delle caratteristiche genomiche. Per ciascuna funzione viene rappresentato il numero di DMR 
sovrapposti rispetto al numero totale di finestre testate per la rispettiva funzione. Rispetto a 
tutte le regioni genomiche, le regioni ipermetilate sono arricchite per satelliti, ripetizioni 
semplici e CGI distali ai promotori senza prove normative note e le regioni ipometilate sono 


arricchite per ripetizioni semplici. 
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Il profilo della metilazione dell'intero genoma mediante MeDIP-seq ha rivelato 66.178 
regioni differenzialmente metilate (DMR; valore q <0,05) in MT cardiaci. Queste 
alterazioni hanno interessato I'1,1% del genoma coperto. In linea con i cambiamenti 
trascrizionali osservati di scrittori e cancellatori della metilazione del DNA, il 71% dei 
DMR corrispondeva al guadagno di metilazione (46.984 regioni ipermetilate contro 
19.194 regioni ipometilate). Al contrario, nei MT epatici, nessun DMR ha superato la 
correzione per test multipli. Inoltre, il grafico PCA (Fig. 2d ) ha indicato una differenza 
tra i gruppi di trattamento nei MT cardiaci ma non nei MT epatici. Insieme, questo ha 
illustrato l'impatto specifico del tessuto del DMSO sul metiloma nella maturazione dei 


MT cardiaci, mentre i MT epatici maturi non sono apparsi inalterati. 


Al fine di analizzare gli effetti regolatori dei DMR, le regioni sono state annotate con 
caratteristiche normative note, come promotori, isole CpG (CGI), siti di legame del 
fattore di trascrizione (TFBS) e diverse classi di elementi ripetitivi (Fig. 4b ). Sia le 
regioni ipo che quelle ipermetilate hanno arricchito classi di ripetizioni specifiche 
come i satelliti (odds ratio di iper/ipometilazione: 4,4/2,4) e le ripetizioni semplici 
(iper/ipometilazione 6,0/9,4 odds ratio), mentre gli elementi ripetitivi in generale non 
sono stati arricchiti. L'effetto del DMSO sulla sovraregolazione dei DNMT e la 
successiva ipermetilazione delle sequenze ripetute era in linea con i risultati 
precedenti nei corpi embrionali di topo '?. Come previsto, anche le regioni 
ipermetilate erano altamente arricchite in CGI senza prove normative, ad esempio non 
sovrapposte con TFBS e distali ai siti di inizio della trascrizione. L'odds ratio di queste 
regioni sul genoma era 5,5. Mentre le CGI in generale avevano un odds ratio di 3,2, 
non c'era arricchimento di CGI che si sovrapponevano a TFBS noti presso i promotori, 
le regioni con il potenziale normativo meglio studiato. Proprio come le regioni del 
promotore, altre caratteristiche come esoni, introni e TFBS sono state influenzate 
dall'ipermetilazione all'incirca alle quote del background genomico catturato (Fig. 
4b). Sebbene abbiamo trovato esempi specifici in cui l'ipermetilazione del promotore 
si verifica in concomitanza con la repressione trascrizionale, non é stata osservata 
alcuna correlazione globale tra la metilazione differenziale del promotore e 
l'espressione genica. Solo 148 dei 1.436 geni (10,396) con promotori ipermetilati erano 
sottoregolati nei MT cardiaci trattati con DMSO. 86 dei geni (696) erano sovraespressi, 


contrariamente all'effetto repressivo previsto della metilazione. 
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Presi insieme, i nostri dati implicavano un'influenza fondamentale del DMSO sulla 
metilazione del DNA nei MT cardiaci, ma non nei MT epatici. Le alterazioni indotte 
hanno interessato preferenzialmente i satelliti e le ripetizioni semplici, nonché i CGI 
distali rispetto ai promotori senza prove normative note. L'espressione dei geni vicini 
non era significativamente correlata alle differenze di metilazione. Concludiamo 
quindi che la repressione trascrizionale da parte dell'ipermetilazione del promotore 
non era l'effetto regolatorio primario dei segni di metilazione alterati, ma i risultati 


indicano un'interruzione globale dei meccanismi di metilazione del DNA. 


Discussione 


Mentre alte dosi di DMSO sono state inizialmente studiate per applicazioni mediche, 
dosi più basse sono ora comunemente usate come solventi per la tossicologia o per la 
crioconservazione di colture cellulari, ovociti ed embrioni. Nel frattempo, gli effetti 
biologici del DMSO sono stati dimenticati o considerati trascurabili. Con l'evoluzione 
della scienza biomedica verso tecniche ad alto rendimento più sensibili e verso nuove 
aree di ricerca, comprese le modifiche dell'epigenoma e il silenziamento genico 
mediato da microRNA, abbiamo valutato gli effetti biologici dell'esposizione al DMSO 


sul trascrittoma e sull'epigenoma. 


La prima panoramica degli effetti cross-omici indotti da DMSO, osservati attraverso 
quantità di entità differenzialmente modificate e analisi PCA, ha indicato una 
maggiore suscettibilità agli effetti DMSO nei MT cardiaci. Le differenze tra i tipi di 
tessuto studiati possono riguardare il fatto che il modello cardiaco contiene 
cardiomiociti umani derivati da iPSC che stanno ancora maturando, mentre il modello 
epatico contiene epatociti umani maturi. Attualmente, non è ancora possibile che le 


20,21 e |a rassegna concisa sui 


cellule iPSC raggiungano un fenotipo simile all'adulto 
cardiomiociti umani derivati da iPSC di Robertson et al . **descrive le iPSC per 
assomigliare meglio ai cardiomiociti fetali (per quanto riguarda la struttura, la 
proliferazione, il metabolismo e l'elettrofisiologia) piuttosto che i cardiomiociti adulti. 
Le nostre osservazioni sugli effetti del DMSO, in particolare le modificazioni 


epigenetiche estreme nelle cellule cardiache simili a quelle fetali, sono particolarmente 
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rilevanti per la crioconservazione durante la tecnologia di riproduzione assistita 
(embrioni, ovociti e spermatozoi). Un'altra differenza da considerare tra i nostri due 
modelli cellulari è il numero di donatori. Infatti, mentre i MT cardiaci sono composti 
da cardiomiociti derivati da IPSC da un singolo donatore (co-coltivati con fibroblasti 
anche da 1 donatore), i MT epatici contengono epatociti maturi da 10 donatori (co- 
coltivati con cellule di Kupffer di 1 donatore). Mentre un modello a donatore singolo è 
più paragonabile a uno in vivosituazione quindi un modello multi-donatore, i MT 
cardiaci a donatore singolo potrebbero essere più suscettibili a bias a causa della 
sensibilità del genotipo rispetto ai nostri MT epatici multi-donatore che hanno la 


capacità di uniformare questi effetti. 


Gli effetti del DMSO sui processi cellulari sono stati analizzati utilizzando i dati 
completi del trascrittoma. Sono stati identificati 2.051 DEG per MT cardiaci e 2711 
DEG per MT epatici. In entrambi i tipi di tessuto, si riscontrano sottoregolazioni per il 
60% dei DEG. L'analisi della sovrarappresentazione del percorso mediante DEG, 
evidenzia effetti DMSO simili in entrambi i tipi di tessuto. Ciò indica azioni robuste di 
DMSO che potrebbero essere estrapolate ad altri tessuti. | cambiamenti osservati nei 
processi cellulari includono cambiamenti nelle vie mitocondriali legate alla produzione 
di ROS e alla generazione di ATP cellulare. DMSO è un noto scavenger di radicali 3 . 
Sebbene questa proprietà del DMSO possa fornire protezione contro alti livelli di ROS, 
a livelli normali o ridotti di ROS, puó ostacolare il metabolismo delle cellule basali 
assorbendo gli elettroni necessari per la produzione di ATP ?^. La conseguente 
diminuzione del contenuto di ATP é stata osservata nei nostri due modelli cellulari 
trattati con DMSO (Figura 2 supplementare ). | MT cardiaci hanno mostrato una forte 
diminuzione nel tempo (diminuzione dell'87% dopo 2 settimane di esposizione), 
mentre il contenuto di ATP dei MT epatici mostra inizialmente una piccola 
diminuzione (diminuzione del 3696 nelle prime 72 ore), dopodiché il contenuto di ATP 
é lentamente recuperato al valore basale livello. Diversi effetti a valle sono rilevanti a 
seconda dell'applicazione biologica del DMSO. Ad esempio, per i saggi cellulari e la 
ricerca tossicologica, i cambiamenti nella quantità di ROS o di energia libera nella 
cellula influiscono sulla capacità della cellula di affrontare gli stress indotti, portando 
cosi potenzialmente a un'interpretazione errata dei risultati in vitrosaggi. All'interno 


della tecnologia di riproduzione assistita, il contenuto di ATP é un buon predittore 
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della vitalità dell'embrione. La diminuzione dell'ATP indotta da DMSO, specialmente 
nella fase di scissione, può indurre effetti a valle che possono interrompere la funzione 
cellulare, la capacità di impianto e lo sviluppo fetale ‘> 2° . La riduzione indotta da 
DMSO dei tassi di impianto e le perdite di gravidanza erano già state osservate da 


modelli animali 2’ . 


Oltre a sovrapporsi nei processi interessati dal DMSO tra i tipi di tessuto, l'analisi del 
percorso ha anche rivelato differenze tessuto-specifiche. La maggiore quantità di 
processi interessati rilevati nei MT cardiaci può indicare una maggiore suscettibilità 
dei MT cardiaci agli effetti del DMSO, sebbene potrebbe esserci anche un pregiudizio 
statistico dovuto alla diversa quantità di DEG utilizzati nell'analisi del percorso. La 


differenza più significativa è il rilevamento dell "organizzazione della cromatina" nei 
MT cardiaci, mentre questo cluster non è osservato nei MT epatici. | MT cardiaci 
mostrano anche più effetti in altri processi che regolano l'espressione genica. | risultati 
dell'analisi dei miRNA hanno rivelato la completa deregolamentazione della biogenesi 
dei miRNA cardiaci e del contenuto di miRNA. Sebbene i cambiamenti nella biogenesi 
dei miRNA epatici fossero minimi, ancora il 18% dei miRNA rilevati è stato 
significativamente modificato dall'esposizione al DMSO. Gli effetti di metilazione 
dell'intero genoma sul trattamento con DMSO mostrano un'ampia variazione tessuto- 
specifica. Sebbene non vi siano differenze significative rispetto ai MT epatici, si 
osservano forti effetti con i MT cardiaci, in particolare rispetto a sequenze ripetitive 
specifiche come satelliti e ripetizioni semplici. Gli elementi ripetitivi del DNA 
rappresentano il 50% del genoma umano e la metilazione in queste regioni sopprime 
la loro mobilità e mantiene la stabilità genomica. | cambiamenti nello stato di 
metilazione delle ripetizioni dopo il trattamento con DMSO potrebbero destabilizzare 
il genoma, implicando conseguenze fenotipiche e patologiche. Infatti, recentemente, 
sono state riscontrate variazioni di metilazione di sequenze ripetitive, in particolare 
elementi di ripetizione satellite, nel tessuto cardiaco cardiomiopatico rispetto al 
normale tessuto cardiaco Sebbene non vi siano differenze significative rispetto ai MT 
epatici, si osservano forti effetti con i MT cardiaci, in particolare rispetto a sequenze 
ripetitive specifiche come satelliti e ripetizioni semplici. Gli elementi ripetitivi del DNA 
rappresentano il 50% del genoma umano e la metilazione in queste regioni sopprime 


la loro mobilità e mantiene la stabilità genomica. | cambiamenti nello stato di 
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metilazione delle ripetizioni dopo il trattamento con DMSO potrebbero destabilizzare 
il genoma, implicando conseguenze fenotipiche e patologiche. Infatti, recentemente, 
sono state riscontrate variazioni di metilazione di sequenze ripetitive, in particolare 
elementi di ripetizione satellite, nel tessuto cardiaco cardiomiopatico rispetto al 
normale tessuto cardiaco Sebbene non vi siano differenze significative rispetto ai MT 
epatici, si osservano forti effetti con i MT cardiaci, in particolare rispetto a sequenze 
ripetitive specifiche come satelliti e ripetizioni semplici. Gli elementi ripetitivi del DNA 
rappresentano il 50% del genoma umano e la metilazione in queste regioni sopprime 
la loro mobilità e mantiene la stabilità genomica. | cambiamenti nello stato di 
metilazione delle ripetizioni dopo il trattamento con DMSO potrebbero destabilizzare 
il genoma, implicando conseguenze fenotipiche e patologiche. Infatti, recentemente, 
sono state riscontrate variazioni di metilazione di sequenze ripetitive, in particolare 
elementi di ripetizione satellite, nel tessuto cardiaco cardiomiopatico rispetto al 
normale tessuto cardiaco in particolare rispetto a specifiche sequenze ripetitive come 
satelliti e ripetizioni semplici. Gli elementi ripetitivi del DNA rappresentano il 50% del 
genoma umano e la metilazione in queste regioni sopprime la loro mobilità e 
mantiene la stabilità genomica. | cambiamenti nello stato di metilazione delle 
ripetizioni dopo il trattamento con DMSO potrebbero destabilizzare il genoma, 
implicando conseguenze fenotipiche e patologiche. Infatti, recentemente, sono state 
riscontrate variazioni di metilazione di sequenze ripetitive, in particolare elementi di 
ripetizione satellite, nel tessuto cardiaco cardiomiopatico rispetto al normale tessuto 
cardiaco in particolare rispetto a specifiche sequenze ripetitive come satelliti e 
ripetizioni semplici. Gli elementi ripetitivi del DNA rappresentano il 50% del genoma 
umano e la metilazione in queste regioni sopprime la loro mobilità e mantiene la 
stabilità genomica. | cambiamenti nello stato di metilazione delle ripetizioni dopo il 
trattamento con DMSO potrebbero destabilizzare il genoma, implicando conseguenze 
fenotipiche e patologiche. Infatti, recentemente, sono state riscontrate variazioni di 
metilazione di sequenze ripetitive, in particolare elementi di ripetizione satellite, nel 
tessuto cardiaco cardiomiopatico rispetto al normale tessuto cardiaco | cambiamenti 
nello stato di metilazione delle ripetizioni dopo il trattamento con DMSO potrebbero 
destabilizzare il genoma, implicando conseguenze fenotipiche e patologiche. Infatti, 
recentemente, sono state riscontrate variazioni di metilazione di sequenze ripetitive, in 


particolare elementi di ripetizione satellite, nel tessuto cardiaco cardiomiopatico 


https://www.nature.com/articles/s41598-019-40660-0 18/39 


24/03/22, 06:09 DMSO induce cambiamenti drastici nei processi cellulari umani e nel panorama epigenetico in vitro | Rapporti scientifici 
rispetto al normale tessuto cardiaco | cambiamenti nello stato di metilazione delle 
ripetizioni dopo il trattamento con DMSO potrebbero destabilizzare il genoma, 
implicando conseguenze fenotipiche e patologiche. Infatti, recentemente, sono state 
riscontrate variazioni di metilazione di sequenze ripetitive, in particolare elementi di 
ripetizione satellite, nel tessuto cardiaco cardiomiopatico rispetto al normale tessuto 
cardiaco?? . Inoltre, è stato dimostrato che il DMSO modifica i profili di metilazione del 
DNA dell'intero genoma in base a specifici loci genici. Ciò potrebbe essere dovuto alla 
drastica stimolazione di DNMTT e DNMT3 . Abbiamo osservato la sovraregolazione 
dei DNMT , ma la letteratura ha anche riportato un aumento dei livelli proteici e 
dell'attività catalitica dall'interazione del DMSO con i loro substrati enzimatici (DNA e 
AdoMet) 19: 29.30. Inoltre, è anche possibile che il DMSO agisca come un donatore di 
metile per indurre l'ipermetilazione ?!. Tuttavia, va tenuto presente che esiste 
un'interazione tra la metilazione del DNA e la configurazione della cromatina, che è 
regolata da più processi rispetto ai soli DNMT. Nel complesso, ipotizziamo che gli 
effetti del DMSO sui cambiamenti della metilazione siano una deregolamentazione 
globale caratterizzata da un'ipermetilazione dell'intero genoma. | meccanismi cellulari 
incaricati di rimuovere le metilazioni dannose in regioni importanti, come TFBS, 
potrebbero quindi aver rimosso i cambiamenti di metilazione con un esito negativo 
sulla sopravvivenza cellulare, lasciando una sovrarappresentazione del DMR negli 
elementi ripetitivi. Questo è in analogia con i tassi di mutazione, che sono più elevati 
per gli elementi ripetitivi poiché sono meno soggetti alla selezione naturale e ai 
meccanismi di riparazione del DNA ?*. Riteniamo che la differenza di metilazione tra 
cellule epatiche mature e cellule cardiache in maturazione possa essere dannosa per 
le cellule crioconservate e in particolare per ovociti ed embrioni. Gli effetti DMSO 
precedentemente riportati su ovociti ed embrioni sono limitati ai cambiamenti 
nell'espressione genica nei modelli animali 27:33 34.35.36, Abbiamo ora mostrato, 
per la prima volta, l'influenza del DMSO sull'epigenoma delle cellule umane in vitro . 
Non è chiaro se le cellule esposte possano riprendersi da esposizioni temporanee a 
DMSO. È stata documentata una maggiore sopravvivenza per gli embrioni dopo il 
ritiro dal DMSO al 2% per meno di 24 ore, sebbene questo non fosse il caso con 
esposizioni più lunghe di 48 ore *’. Sebbene i cambiamenti di metilazione abbiano 
una natura adattiva, possono anche essere persistenti. Dalla ricerca condotta sulla 


coorte olandese della carestia, è noto che il periodo periconcezionale (prima del 
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concepimento fino all'inizio della gravidanza) è cruciale per stabilire e mantenere 
segni epigenetici. Durante questo periodo i segni epigenetici presenti sul DNA 
parentale vengono rimossi per produrre uno zigote totipotente che viene 
riprogrammato dopo l'impianto 3”. Le esposizioni durante questo periodo possono 
indurre differenze epigenetiche persistenti, che possono provocare malattie più tardi 
nella vita o eventualmente essere ereditate dalla prole per indurre effetti 
transgenerazionali °°: 79: 4°: 41,42. È già noto che i bambini concepiti attraverso le 
tecnologie di riproduzione assistita tendono a nascere pretermine, hanno un basso 
peso alla nascita e hanno anche un rischio significativamente maggiore di difetti alla 


43,44. Inoltre, alcuni esiti di salute a breve termine sono stati leggermente 


nascita 
inferiori nei bambini concepiti con la fecondazione in vitro, sebbene nel complesso i 
risultati siano stati positivi 42 4° , Per quanto riguarda la salute e lo sviluppo a lungo 
termine, sebbene i dati disponibili siano limitati, sono stati indicati i fattori di rischio 
cardiovascolare e metabolico che alla fine possono sfociare in malattie 
cardiometaboliche croniche ^^. Inoltre, sono stati dibattuti anche i rischi epigenetici, 
concentrandosi principalmente sulla metilazione e sull'imprinting. Per i disturbi 
dell'imprinting più studiati, la sindrome di Beck-with-Wiedemann e la sindrome di 
Angelman, l'incidenza era maggiore per i bambini nati dopo le tecnologie di 
riproduzione assistita rispetto al concepimento naturale ^! ‘46. 47. 48,43,50 Sebbene 
ci siano molti fattori che contribuiscono a queste anomalie (diminuzione dell'infertilità, 
età più elevata dei genitori, interferenza tecnica con un processo biologico, ecc.) Gli 


effetti indotti dal DMSO possono essere uno di questi. 


Nell'ultimo decennio, i miglioramenti apportati alla tecnologia riproduttiva artificiale 
hanno notevolmente aumentato il tasso di successo dei trattamenti di fecondazione in 
vitro. In particolare, le migliori tecniche di crioconservazione (transizione alla 
vetrificazione) hanno aiutato notevolmente. Sebbene non sia stato ancora dimostrato, 
vengono riportati tassi di successo aumentati quando si impiantano embrioni 
precedentemente congelati ^! 52: >3 , Ciò ha indotto la tendenza del "congela tutto" 
nella clinica, che aumenta notevolmente la necessità di comprendere le conseguenze 
della crioconservazione a livello molecolare. Inoltre, mentre la conservazione degli 
ovociti congelati era inizialmente destinata alle donne minacciate dalla diminuzione 


della fertilità dovuta a malattie, nella nostra società incentrata sui portatori, è 
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diventata anche un metodo per prevenire l'infertilità legata all'età per le donne che 
vogliono ritardare la maternità ?^. Lo sviluppo più recente nella tecnologia di 
riproduzione assistita è l'emergere di preoccupazioni sugli effetti DMSO. Sebbene ci 
siano ancora controversie su questo argomento, il campo si sta muovendo verso 
metodi privi di DMSO. Tuttavia, privo di DMSO non significa automaticamente sicuro 


a causa dell'uso di altri crioprotettori ??: 36.37, 


Conclusioni 


Il nostro studio mette in evidenza la capacità del DMSO di indurre cambiamenti nei 
processi cellulari sia nelle cellule cardiache che epatiche, ma più gravemente, indurre 
alterazioni nel miRNA e nel panorama epigenetico nel modello cardiaco in 
maturazione 3D. | cambiamenti nei processi cellulari possono avere conseguenze per 
le conclusioni tratte dai saggi cellulari e quindi anche in qualsiasi applicazione di 
questi risultati (ad es. conclusioni false negative sulla tossicità dei farmaci). Inoltre, i 
cambiamenti estremi nei miRNA e le alterazioni del panorama epigenetico possono 
rappresentare una minaccia, soprattutto per la tecnologia di riproduzione assistita. 
L'ipermetilazione dell'intero genoma indotta dalla deregolamentazione globale dei 
meccanismi di metilazione, specialmente quando colpisce geni importanti nello 
sviluppo, può avere conseguenze negative direttamente, più avanti nella vita o forse 


nelle generazioni successive. 


Nel complesso, l'uso di DMSO dovrebbe essere evitato ove possibile. Tuttavia, per il 
momento, il DMSO è indispensabile nelle applicazioni biotecnologiche. In questi casi, 
dovrebbero essere considerati gli effetti che il DMSO può avere e la concentrazione 
dovrebbe essere mantenuta il più bassa possibile, perché anche a basse 


concentrazioni il DMSO non è inerte. 


Metodi 
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Campioni 
In questo studio sono stati utilizzati microtessuti cardiaci umani 3D InSight ™ 
(InSphero, per beta-tester), contenenti una co-coltura di circa 4000 cardiomiociti 
umani derivati da iPSC da una donatrice caucasica senza fenotipo di malattia noto e 
1000 fibroblasti cardiaci da un Maschio caucasico di 18 anni. Le MT sono state 
coltivate in 3D Insight "M Human Cardiac Microtissue Maintenance Medium 
(InSphero). Inoltre, sono stati utilizzati 3D InSight '™ Human Liver Microtissues 
(InSphero), circa 1000 epatociti umani primari (pool multi-donatore di 5 maschi e 5 
femmine di età compresa tra 7 e 59 anni) in co-coltura con 1000 cellule di Kupffer 
umane primarie da un caucasico 27 anni di sesso non dichiarato. Le MT sono state 
coltivate in 3D Insight "MMezzo di mantenimento del microtessuto del fegato umano 
- AF (InSphero). 


Esposizione 

| microtessuti (MT) sono stati esposti solo al mezzo (= non trattato) o al mezzo 
contenente lo 0,1% di DMSO per due settimane. Il mezzo è stato cambiato tre volte al 
giorno, in accordo con il disegno sperimentale della nostra attuale ricerca °° . Durante 
l'esposizione di 2 settimane, ci sono stati un totale di 7 punti temporali di 
campionamento a 2, 8, 24, 72, 168, 240 e 336 ore, che sono stati campionati in 
triplicato (totale 21 campioni per condizione). Per campione, sono stati utilizzati 36 MT 
cardiaci o 54 MT epatici per ottenere quantità sufficienti di RNA per il 


sequenziamento. 


Preparazione del campione di sequenziamento dell'RNA 

L'RNA totale è stato isolato utilizzando il kit universale AllPrep DNA/RNA/miRNA 
(Qiagen). Il campione è stato impoverito di RNA ribosomiale utilizzando il kit di 
rimozione dell'rRNA Ribo-Zero Gold (umano/topo/ratto) (Illumina ®) e preparato per 
il sequenziamento con il kit di preparazione della libreria di RNA totale SENSE 
(Lexogen). Dopo la preparazione della libreria, la qualità è stata valutata su 
TapeStation Agilent 4200 e la concentrazione della libreria è stata determinata da 
Qubit ™ prima del sequenziamento (paired-end 100 bp) su HiSeq2000. Sono stati 


ottenuti due file per campione, divisi in letture sinistra e destra e le sequenze 
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dell'adattatore Lexogen (le prime 12 basi di tutte le letture) sono state rimosse 
utilizzando Trimmomatic ?? (v.0.33). La qualità dei dati di sequenziamento è stata 
verificata utilizzando FastQC ©°(v.0.10.1) prima e dopo la rifilatura. Per i campioni 
cardiaci, nessun campione è stato scartato a causa della scarsa qualità. Due campioni 
sono stati scartati a causa del basso numero di letture (UNTR_002_3 e UNTR_240_3) e 
l'analisi dei componenti principali (Figura 1 supplementare ) ha rivelato 4 valori 
anomali che sono stati rimossi (UNTR 008 3, UNTR 168 2, UNTR 2402 e 
DMSO 336 2). Infine, poiché é rimasta solo una replica di 240 ore di campioni non 
trattati, tutti i campioni di questo punto temporale sono stati rimossi dall'analisi 


statistica. Per i campioni epatici, nessun campione é stato scartato. 


Analisi dei dati dell'RNA 

| conteggi delle letture per i geni sono stati ottenuti allineando le letture al genoma di 
riferimento (Genome Reference Consortium Human Build 38, GRCh38.p7) utilizzando 
RSEM ©! (v.1.2.28) con il paired-end e Bowtie2 © (v.2.2.6 ) opzione. Successivamente, i 
conteggi letti per ciascun gene sono stati utilizzati per la determinazione di geni 
differenzialmente espressi utilizzando il pacchetto DESeq2 9? R (v.1.16.1). Qui, il 
design é stato impostato in base all'esposizione (UNTR o DMSO), le altre impostazioni 
sono state mantenute sui parametri predefiniti, ad eccezione del conteggio minimo, 
che é stato impostato su un conteggio medio di lettura di uno su tutti i campioni. 
Infine, i geni ribosomiali rimanenti sono stati filtrati dal set di dati per garantire il 


completo esaurimento ribosomiale dei dati. 


Preparazione del campione di sequenziamento del MiRNA 

Un'aliquota dell'RNA totale isolato é stata selezionata in base alla dimensione e ligata 
utilizzando il TruSeq Small RNA Library Prep Kit (Illumina 9). Dopo la preparazione 
della libreria, i campioni sono stati sequenziati (paired-end 100 bp) su HiSeq2000. Le 
sequenze dell'adattatore sono state tagliate di 9^ e le letture comprese tra 16 e 35 bp 
sono state mantenute per l'analisi. Tre campioni sono stati scartati a causa del basso 
numero di letture (cardiaco: UNTR, 168 1, UNTR 168 2 ed epatico: DMSO 002 3). 
Inoltre, l'analisi delle componenti principali (Figura 1 supplementare ) non ha rivelato 


valori anomali. 
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Analisi dei dati MIRNA 
L'analisi dei dati di MIRNA è stata eseguita in modo simile come descritto in 
precedenza °°, ad eccezione dell'uso di mirDeep2, che non ottiene informazioni 
aggiuntive quando si utilizzano campioni umani. In breve, PatMaN °° (v.1.2.2) è stato 
utilizzato per allineare le letture tagliate al genoma umano senza consentire 
discrepanze o lacune. Per ottenere conteggi di lettura completi per le specie di miRNA 


3' e 5', l'output della mappatura è stato analizzato. 


proteomica 


Anche i dati di proteomica sono stati ottenuti e analizzati. | metodi ei risultati sono 


inclusi nei dati supplementari. 


Sequenziamento MeDIP 


Per la preparazione delle librerie MeDIP-Seq, il protocollo MeDIP 97 


a basso input è 
stato modificato. Il DNA è stato frammentato a 100-200 bp utilizzando il sistema 
Covaris S2. A causa della minore resa del DNA per i campioni cardiaci di DMSO allo 
0,1%, i triplicati sono stati raggruppati prima della frammentazione. La riparazione 
finale e l'accodamento sono stati eseguiti utilizzando il kit di preparazione della 
libreria NEBNext® Ultra ™ per Illumina® (NEB), gli adattatori sono stati legati con 
NEBNext® Ultra MII modulo di ligazione (NEB) e i campioni sono stati purificati 
utilizzando perline Agencourt® AMPure® XP (Beckman Coulter). | frammenti metilati 
sono stati catturati utilizzando il kit MagMeDIP (Diagenode). In breve, il DNA 
denaturato è stato miscelato con l'anticorpo anti-5-meC e catturato utilizzando sfere 
magnetiche. L'efficienza di acquisizione è stata determinata mediante qPCR contro 
frammenti di Lambda-DNA addizionati in campioni di libreria di precattura e post- 
acquisizione. Le librerie sono state amplificate in una fase finale della PCR utilizzando 
primer TruSeq con codice a barre. La qualità è stata valutata su Agilent Bioanalyzer 


2100 e la concentrazione della libreria è stata determinata mediante Qubit ™ e qPCR. 


Analisi dei dati MeDIP 


Al fine di ottenere una copertura completa dell'intero genoma, i campioni triplicati 


che sono stati sequenziati individualmente sono stati uniti prima dell'allineamento. Le 
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letture del sequenziamento MeDIP sono state allineate al genoma di riferimento 
GRCh38 utilizzando bwa versione 0.7.15-r1140 88 e analizzate in 250 finestre di base 
utilizzando il pacchetto R/bioconduttore QSEA €? (v.1.4.0) con parametri standard. 
All'interno di QSEA, l'arricchimento MeDIP è stato calibrato con misurazioni dell'array 
di metilazione 450k di epatociti primari (GSM999339) e miociti cardiaci (HCM, 
GSM999381) da ENCODE ”°, rispettivamente per i micro tessuti epatici e cardiaci. A 
tal fine sono stati calcolati i valori beta dei campioni di calibrazione con il pacchetto 
R/bioconduttore Minfi ’'(v.1.24.0), le posizioni genomiche delle sonde dell'array sono 
state mappate da GRCh37 a GRCh38 utilizzando lo strumento a riga di comando 
UCSC liftOver "^ e si è calcolata la media delle sonde all'interno di 250 finestre di base. 
Le regioni metilate differenzialmente ottenute da QSEA sono state annotate con 
informazioni su gene, esone e promotore (sito di inizio della trascrizione +/- 2 
kilobasi) da RefSeq, ENCODE TFBS e isole CpG basate su modello, tutte ottenute 
tramite il browser della tabella UCSC. Poiché ENCODE TFBS non era disponibile per 
GRCh38, le posizioni genomiche sono state mappate da GRCh37 utilizzando lo 


strumento liftOver. 


Disponibilità dei dati 
| set di dati generati e analizzati durante il presente studio sono disponibili nel 
database BioStudies ( http://www.ebi.ac.uk/biostudies ) con i numeri di accesso S- 


HECA1, S-HECA5, S-HECA9, S-HECA 18, S -HECA33, S-HECA34, S-HECA47, S-HECA139, 
S-HECA158, S-HECA361, S-HECA401, S-HECA403. 
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